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HIGHER-ORDER NUMERICAL SOLUTIONS USING CUBIC SPLINES 


S. G. Rubin and p 0 K„ Khosla 

Polytechnic Institute of New York 
Farmingdalc, New York 

ABSTRACT 

A cubic spline collocation procedure has Recently been 
developed for the numerical solution of partial differential 
equations o In the present paper, this spline procedure is 
reformulated so that the accuracy of the second-derivative 
approximation is improved and parallels that previously ob- 
tained for lower derivative terms. The final result is a 
numerical procedure having overall third-order accuracy for 
a non-uniform mesh and overall fourth-order accuracy for a 
uniform mesh. Solutions using both spline procedures, as 
well as three-point finite difference methods, will be pre- 
sented for several model problems. 

I „ INTRODUCTION 

. 1 2 

In a recent study Rubrn and Graves ' have presented a 
cubic spline 3, ^ collocation procedure for the numerical solu- 
tion of partial differential equations. This technique ex- 
hibits the following desirable features: (1) The governing 

matrix system is always tridiagonal so that well-developed, 
and highly efficient inversion algorithms are applicable; 

(2) cubic spline interpolation leads to second order accuracy 
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for second derivatives, e.g., diffusion terms in the Navier- 
Stokes equations. This order of accuracy is maintained oven 
with rather large non-uniformities in mesh width? (3) first 
derivatives or convection effects are fourth-order accurate 
for a uniform mesh and third-order with mesh non-uniformity? 

(4) dorivseive boundary conditions can in many cases be applied 
more accurately and with less difficulty than with conven- 
tional finite-difference schemes? (5) a simple two-point 

* i 

relationship exists between the spline approximation for the 
first and second derivatives? and (6) unlike finite-element 
or other Galerkin (integral) methods, which are generally not 
tridiagonal, the evaluation of Sarge numbers of quadratures is 
unnecessary. 

1 2 

Solutions have been obtained for a number of problems ' 
with explicit,, implicit and spline alternating direction 
implicit (SADI) temporal or spatial marching procedures. 

Moreover, for the viscous and potential flow problems consid- 
ered, it was found that with the spline procedure there was 
no particular advantage gained with the equations in divergence 
form. In some recent studies it has been found that the 
divergence form may be desirable with flux boundary conditions. 
These results are described later in this paper. 

Agreement of the spline solutions with exact analytic 

results and very accurate finite-difference solutions obtained 

12 12 

with a very fine mesh has been quite good ' . All comparisons ' 
with conventional three-point finite difference formulations 
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demonstrate the improved spline accuracy associated with 
(i) the higher-order convection approximation, (ii) the treat- 
ment of derivative boundary conditions, or (iii) the higher - 
order accuracy of spline second derivatives (diffusion) when 
specifying a non-uniform mesh. Solutions for the Burgers 
equation, the two-dimensional diffusion equation and the in- 
compressible viscous flow in a driven cavity are found in 
Refs. 1 and 2„ 

In the present paper, the cubic spline procedure is 

reformulated so that the accuracy of the second-derivative 

approximation is improved and parallels that obtained for the 

lower derivative terms „ The final result is a combined spline- 

finite difference numerical procedure having overall third- 

order spatial accuracy for non-uniform meshes and overall 
* 

fourth-order spatial accuracy with a uniform mesh. In order 

to differentiate the two spline procedures, we shall designate 

1 2 

the original spline formulation ' as spline 2 and the improved 
formulation presented here as spline 4. 

hs shown in sections II and III, the cubic spline colloca- 
tion procedure involves a third-order interpolation polynomial 
with the function and the second (or first) derivative of the 
function as unknowns at each mesh point. Continuity of the 
first (or second) derivative leads to the tridiagonal system 
of equations to be considered. In section IV, it is shown how 
the familiar central difference second-order accurate finite- 
difference theory results from a quadratic spline interpolation 
procedure. Using the earlier spline designation, the finite- 
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difference theory is classified as spline 1. 

Recently, several higher-order finite-difference schemes 
with similar properties have been proposed; i.e,, the functions 
and derivatives are considered unknown at each mesh point, or 

the functions are collocated at three points instead of one. 

. 5 6 

The methods which have been termed Hcrmitian finite-difference ' , 

/ 7 8 9 

Padc approximation ' or compact differencing , and Mehrstellung 

have been developed for a uniform mesh and have somewhat lower 

truncation errors than the five-point pentadiagonal fourth-order 

finite difference procedure . As with the spline formulation, 

they remain of tridiagonal form* 

In a recent study the authors' 1-0 have examined these pro- 
cedures, as well as a fourth-order spline-on-spline method, and 
found them to be, in fact, identi cal; i.e. , any one can be 
derived from any of the others. Moreover, these procedures 
have been reformulated' 1 ' 0 so as to apply to non-uniform mesh 
systems as well. As with spline 4, these finite-difference or 
spline-on-spline methods are fourth-order with a uniform mesh 
and third-order with a non-uniform mesh. The relative advantage 
and/or disadvantages of these procedures, over spline 2 and 
spline 4, are discussed in Ref. 10. The main differences are 
handling of the boundary conditions, the relationship between 
the approximations for the convection and diffusion terms, and 
the truncation errors. The truncation errors for first deriv- 
atives are identical. The truncation error for the second- 
derivative to be discussed later for spline 4 is 50% smaller 
than that found with the higher-order finite-difference or 
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splinc-on~splinc collocation formulae. 

In order to evaluate the spline procedures, the trunca- 
tion errors, stability limitations and effects of boundary 
conditions will be discussed. Spline 2 is reviewed in soction 
IX, spline 4 is introduced and discussed in section III, and 
spline 1 is presented in section IV, The stability conditions 
for all methods are outlined in section V, Solutions using 
both spline procedures, as well as a three-point f inite-dif forence 
method, are presented for several model problems in section VI. 
Both uniform and non-uniform meshes arc considered. In each 
case the analytic solution or a very accurate numerical solu- 
tion is available for comparison purposes. The problems to be 
considered include (1) a boundary laycr-like solution of Laplace's 
equation, where a spline relaxation method is applied, (2) poten- 
tial flow over a circular cylinder with a spline successive 
approximation procedure, (3) the weak shock solution for the 
nonlinear Burgers equation by a two-step expD.icit or an implicit 
spline integration, (4) divergence and non-divergence solutions 
for the linear Burgers equation with flux and other derivative 
boundary conditions, (5) the impulsive motion of right angle 
corner (Rayleigh problem) with SADI, (6) t’ae solution of the 
two-point boundary value problem describing similar boundary 
layer behavior, and (7) non-similar constant pressure boundary 
layer solutions for large Reynolds number using physical variables. 
The results will be summarized in section VII „ 
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IX. SPLINE 2 - REVIEW OF CUBIC SPLINE THEORY 


Consider a mesh with nodal points such that 
a =x o < Xl .;x 3 ... <x N <x N . Hlt =h, 


and with 


h i ex i- x i-i 


> 0 . 


Consider a function u(x) such that at tho mcsh ( points x^, 

u(x^)cu^, Tho cubic spline is a function S^(u^, x)=S^(x) which 

is continuous together with its first and second derivatives on 

the interval [a,b], corresponds to a cubic polynomial in each 

sub-interval x. , <x<x., and satisfies S . (u . ?x. ) «=u. , In the 
r-l— — x A 3.' x' x 

usual spline terminology, spline 2 is defined as a cubic spline 
of deficiency one, since all but one of the three polynomial 
derivatives are continuous. 

If u(x) and its derivatives are continuous, it has been 
shown that the spline function S^(x) approximates u(x) at all 
points in [a,b] to fourth order in max li. . First and second 
derivatives of S^(x) approximate u'(x) and u"(x) to third and 
second order, respectively. See Ahlberg, Nilson and Walsh for 
detailed proofs of convergence. 

If S ^ (x) is cubic on [x iM ^,x^], then in general. 


s ;< x > =M i-i ( ¥-> - v- 


x-x 


h i 


x— 1^ 


where S^x^). 

Integrating twice leads to the interpolation formula on 
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s A (x)t vl 


(Xi-x) a _ 

~~6h. + M i — ?iTT — 

1 1 


< u i-r ~s — 1 




M i h ? 

t-- +(v -?-> 




h. 

i 


(la) 


The constants of integration have "been evaluated from s ^( x j_)” u i 

and s ^( x i - x^“ u d“l‘ S A^ X ^ on ^ x i ,x i*l*l^ is obtained with i+1 
replacing i in (la) , 

The unknown derivatives M. are related by enforcing the 

4 * 

continuity condition on S^(x) . With s ^( x f[) =m x on [i-l/i] and 
S'(xJ)=mt on WQ rec ^ uire m i ” n i “ m i* We £or 

1=1 $ m • O # f N | 


h i h i +h i+l 

“ M. - + ~- r— M.+ 
D 1-1 3 1 


i-l-1 


M i+1= 


(u i+l“ u i ) (u i~ u i-l } 


h i+l 


h i 


(lb) 


Additional spline relationships that are easily derived are 
listed below: 


r 1 - m. n +2(r~- + \ ~ — )m. + 
h i 1 - 1 h i h i+l 1 


3 (u 


l+l 


m i+r 


i+l 

h? 

n i+l 


-u . ) 

i' 






(lc) 


m i+r m i = 


h i+l 


(M i +W i+l> 


(Id) 


h i h i 

m i c T M i+ it M i-i + hj 


u . -u . , 

1 1-1 


^i+l b i+i 

m .= - -±: —A M . +i * h 


u . . , -u . 
1+1 1 


i+l 


(la) 

(If) 
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2m. - 4m. u.-u. . 

ii 3*^1 i 3 . 1 

v + ^ r~ 


2m i+i ,, u i+r u i 

h +6 "IT 3 

n i+l i+1 


(lg) 

(lh) 


Eq3 0 (lb) or (lc) lead to a system of N equations for the N+2 
unknowns M. or m . , respectively. The additional two equations 
are obtained from boundary conditions on m Q0 or M q , 

The resulting tridiagonal system for or m^ is diagonally 

3 

dominant and solved by an efficient inversion algorithm . 

1 2 

Spline 2 for Solving Partial Differential Equations ' 

If the values u^ a~e not prescribed but represent the solu- 
tion of a quasi-linear second order partial differential equa- 
tion, u.=f(u, u , u ), then an approximate solution for u. can 
be obtained by considering the solution of 




<Vi 




This formulation is designated spline 2. If the time derivative 
is discretized in a simple finite-difference fashion, we have 

u n+? _ u n 

— ~ = (1-0) f n +0f n+1 ; (2a) 

0=0, explicit; 0=1, implicit; 0=^, Crank-Nicolson 0 For the 
explicit integration the stability limitations are quite severe, 
see Refs, 1,2 and section VI. Therefore a two-scep procedure is 
considered and is given as: 
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Step 1: 


n+1 
x 



At 




* 


Step 2: 


,n+l 

X 



At 


|nH*l 


Example ; 

Consider the linear Burgers equation 

u t+uu x «v u xx ; U=u(x,t) ; v=v(x,t) . 


(2b) 


(3a) 


With (lb) and (lc) we obtain a system of 3N equations for 
3(N+2) unknowns (see Refs. 1/2 for further details on the 

i *J* 

derivation) . The system (2) can be written as 


where 


»i^X’l +X+0 l’lS- 0 l V l«it< I .«l + l3 - (3b) 


R i= 


0 

•1/h. 

3/h| 


0 

0 

i/h, 


Yi n 

h i /6 


B r 


a. 


( 1+1 /a)/h i 
(--3 (1-1/ a s ) /h® 


a. 


2(l+l/a)/h i 


a 3 i 
(o+l)h./3 


c i B 


0 

1 


i+1 


-3 

hT 


0 
0 

1 


Ys 

h,^ 

3 .+ 1 

"6 

0 


L i+l “i+1 


i 

It is possible to treat the viscous terms (M^) implicitly ( G=l) 
and the convection terms explicitly. As shown in Refs. 1,2, the 
stability of the two-step procedure f'or viscous flows is improved. 

number of variations on this system can be derived with the 
relations (1) . 
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and 



L 0 


Pi pa - 

0 0 

0 o . 

0 6i i 

0 0 

0 0 - 


T 

V. = [u., m., M.] , 


(3c) 


(fi=h i+l /h i 7 Yi=Y 3 =6!=0 

a Q =l 7 cti = 0uV 1 1 At ; a 3 = - Ov^ t * 1 At ; 


p o =i ? P1= -d-0) u? At ? p 3 = (1- 9) v^At . 


(3d) 


A significant advantage of the spline 2 formulation is that with 

expressions (1) it is possible to reduce the 3x3 matrix system 

(3) to a scalar set of equations for alone. The details of 

this reduction process are found in Refs 0 1,2. 

For equations with two space dimensions such that u,=f(u,u , 

u y' u xx ' u yy)' a spline alternating direction .implicit (SADI) 

1 2 

procedure has been presented by Rubin and Graves ' . A spline 
successive approximation method can also be simply formulated. 
Both techniques are discussed later in this paper where several 
example problems are presented. * 

Truncation Error 

For interior points,, the spatial accuracy of the spline 
approximation can be directly estimated from the formulas (lb) 
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and (le) or (If). Expanding m. , M. and u. in Taylor scries 

t J* *L X 

and assuming the necessary continuity of derivatives for u(x,y), 
we obtain# with cFh^ + j/h^, 

(u xx ) i =M i + (u iv ) i h| ( 0 -+D /12 ( o+l) 

-(u V ) .h?(ff~l) (2a s +5cf+2)/180 

X X 

-(u vi ) .hKa®/360H-(a-l) 3 (7cr a -2a+7)/1080‘J 

1 X 

+ 0 (hj?) , (4a) 

and 

(u x ) i !=m i +(uiV) i h i 0(a “ 1)/24 + 

t 

+ (u V ),ht a[l+o(oHl) 3/180 + 0(hf ) . (4b) 

XX X 

11 12 

Fyfe and Daniel and Swartz have presented similar relations, 
for constant h^, in their collocation analysis of cubic splines 
for the solution of two point boundary value problems. 

Therefore* the spline approximation with a non-uniform mesh 
is second-order accurate for and third-order for itk . For a 
uniform mesh rru becomes fourth-order with remaining second- 
order accurate „ in the next section a finite-difference expres- 
sion for (u lv )^ is used to increase the accuracy of It and hence 
the overall accuracy of the procedxire. With this modification 
this formulation will be termed spline 4. 

If (lc) is used to evaluate the truncation error for hk, the 
constant 24 in the second expression on the right-hand side 
becomes 72. For the uniform case, (4b) is recovered in all 
cases. 
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IIIo SPLINE 4 - DERIVATION AND DISCUSSION 


In order to improve the overall accuracy of the spline 2 
formulation, it is necessary to reduce the order of the trunca- 
tion error for (u ) . in (4a) , Although a number of procedures 

XX X 

are possible, we have chosen a very simple modification, whereby 

the error term in (4a) for (u iv ).£ is approximated by a three- 

point discretization for This approximation is first-order 

accurate with a non-uniform mesh and second-order with a uniform 

mesh. Therefore the spline approximation for (u ) . is improved, 

and parallels that for (u ).; ice., third-order accuracy is 

achieved for a non-uniform grid and fourth-order accuracy for 

uniform mesh. This improvement leads to what is termed spline 

4, or a quintic spline of deficiency three. 

10 

In a separate study this procedure is described in some- 
what greater detail and is also applied to the Hermitian or 
'5-9 

Pade fxnxte-dxf ference and splme-on-splxne methods in order 
to develop these procedures for non-uniform grids. 

The development of spline 4 is as follows: The expression 

(4a) can be rewritten in the form 

(u xx ) i aM i +h i o(o+1) MM^./12-K) (<*~ 1)h i' h !) ' (5a > 

where A= (1+a 3 ) /<j(l+a) 3 . 

15 

The familiar three-point discretization formula is 
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(5b) 


< M *x>r ^iThT^i+i' (1+a) M i +aM i-i ] 

- ( a~l) h^ { M xxx ) j/3 “h ? ( 1+ a 3 ) (M iv ) ./12 (1+a) +0(hj) . 

Therefore, (4a) or (5a) becomes 

(u xx> r M i + 1 L/<i > ( M i+1- ( 1+ °> M i+ 0M i-x^ 

-7h*? (H-a 3 ) (a-1) (u v ) ./180-M(u vi ) , [a 3 /360 

4« X J. iii • 

+ (c~l) 3 (7a 3 -2a+7)/1080 3 + 0(hf) . (5c) 

• 1 * 

With (4b), 

(u x ) i =m i -i-° (<c~l)h| f ,h*) t (5d) 

and we obtain a uniform higher-order approximation termed 

•h t 

spline 4. When a=l„ spline 4 xs fourth-order accurate and the 

truncation error of (5b) is smaller than that obtained with 

10 * 5-9 

spline~on-spl.ine „ I-lermitxan or Pade methods , which are in 

turn smaller than the error obtained with five-point finite- 

different discretizations. 

In the spline 4 procedure the relations (lb-lh) still apply; 
however, the interpolation polynomial is no longer applicable 
as spline 4 represents a higher-order interpolation. This point 
is discussed in greater detail in Ref. 10. The governing system 

’+ 

It is possible to apply (5b) to (4b) to make (u ) . fourth-order 
even with a non-uniform mesh, see Ref. 10. x 1 

+ 9 

Higher-order procedures, e.g., spline 6, can be derived in a 
similar manner 10 , and spline 2 is recovered from spline 4 with 
A set equal to zero. 
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remains tridiagonal, unlike spline 2, where the system can be 
reduced to that for alone, the appearance of off-diagonal 
terms in (5b) restricts the reduction process to a 2x2 
system in (u^,M^) . 

For the linear Burgers equation the system is still of 
the form (3b) with 

Yi= “ v ? +1 eAtaA/6 ; - v? +1 0AtA/G t 

X 1 

a s" " V^ +1 0At (l- (1+a) A/6) ? 

p 3 = (1-0) Vi +1 At(l-{l+a) A/6) 7 

6i= vJ +1 (l-e)At A/6 . (6) 

All other entries in (3c, 3d) are unchanged. 


IV. FINITE-DIFFERENCE THEORY /SPLINE 1 


If the procedures given previously for spline 2 and spline 4 
are repeated for a quadratic polynomial interpolation with both 
derivatives continuous (a quadratic spline of zero deficiency) , 
we find on 


S A (x)=u i (x-x._ :L )/h+u i _ 1 (x.-x)/h+(u i -u._ 1 -hm i )(x-x i _ 1 ) (Xj-xl/h 3 , 


where 


and 


S A (x i)" u i' V X i-l ,=U i-.l' S A <X i )=n V 


M - » S'(x7)= - 2(u i -u i _ 1 -m i h)/h^ 
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°n Cxi.x^], with S A (* i+l )-n i+1 , S A ( Xi )»u. and S^(x.)=m. 
we obtain 


M 1 - - 


- 2 < u i- u i + i +ra i h >/ h3 


From the continuity of the second-derivative 



and therefore 


m 


L= < u i+l- u i-l ,/2h 


The expression for becomes 


M i =M i = (u i+i -2u i +u i-i )/h ‘ 


Therefore the quadratic spline of zero deficiency leads to the 
central difference expressions. Similar results are obtained 
if a non-uniform mesh is considered .' 1 ' 0 


Vo STABILITY 


For the linear Burgers equation (3), with u,v constant, the 
interior point stability can be assessed with the von Neumann 
Fourier decomposition of the system (3) for h^=h=constant. 

W ith V^ +r =V^expI w ( x ± +r h) , I=(-lp, (3) becomes 


T:v*}' hl ~p. t v n or , 

1 X X X X X 


-*1 

where G.«T. P. is the amplification matrix. The von Neumann 

J* J. i 

condition necessary for the suppression of all error growth 
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I 


requires that the spectral radius p(G^)<l. The eigenvalues of 
are 

For the one-dimensional equation (3), three numerical pro- 
cedures were considered: (1) convection (m^) and diffusion (,M^) 

explicit, (ii) convection explicit, diffusion implicit (two 
steps required for inviscid stability), and (iii) diffusion and 
convection •implicit. With explicit convection, (i) or (ii), 
both divergence and nondivergence forms of the ‘equations have 
been evaluated in Refs. 1 and 2. 

The stability conditions imposed on these schemes is determ- 
ined from 

uj < i 

(i) Explicit convection and diffusion: 0=0 in (2a, 3). 

Spline 2’ Lf2 : | | 3 = (l-6p (1-coscp) (2+coocp) ” 1 ) 3 +c 3 $ 3 <1, where 

P=vAt/h a , c=uAt/h, $=3sincp/ (2+coscp) , cp=uh. Necessary stability 
limits are 

(a) |3 < 1/6 , 

(b) c < (3)"^ , 

(c) R c =c/p=uh/v < 2(3)^ . . (7a) 

These results are more restrictive than the limits found for the 

. . 15 

forward time central space explicit finite-difference method 

or spline 1, which are 

(a) 3 < 1/2, (b) c < 1 , (c) R c < 2 . 
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Spline 4 : |l | a » (l- (5+coscp) (1-coscp) p/U+coscp)) 15 *!* (3csincp/ 

(2+coscp)) a «l» so that necessary stability limits arc 

(a) 3 <1/4, (b) c<( 10 )Vg , (c) R c <(40) i5 /3 . (7b) 

Once again these conditions arc somewhat more restrictive than 

those obtained with second-order finito-dif ferences (7b) . The 

k 

Pade finite-difference limitation c<(6) y6 is even more restric- 
tive, see Appendix of Rcfs„ 1 and 2. It is significant that in 
all cases the explicit method is unconditionally unstable for 
inviscid flow? i«c„, p=0. 

(ii) Two-step explicit integration (2b) : 

This procedure, which alleviates the inviscid instability 
found in (i), is a two-step predictor-corrector method (see 
Refsp 1 and 2) and is similar to the Brailovskaya two-step 
finite-difference technique. For 0=0, We obtain 

c < $^ n =C(2+cosfp) (Ssincp) -1 ] min =(3) . (7c) 

This result is more restrictive than the c<l CFL condition found 

for the Brailovskaya finite-difference method. 

For 0/0, the effect of diffusion when treated implicitly is 

to improve the inviscid stability limitation. For u-0, the 

1 2 

method is unconditionally stable e , Since the convection terms 
are unchanged, spline 4 has the same stability condition. 

(iii) Implicit convection and diffusion: 

The spline 2 and spline 4 procedures are unconditionally 
Stable if 0>l/2 in Eqs. (2), (3). 
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(iv) SADI: 

In Ref, 1, the interior point stability analysis is extended 
to the two-dimensional SADI procedure; unconditional stability 
is demonstrated. 

Although the implicit procedures lead to unconditionally 
stable formulations, as with finite-difference methods, the tri- 
diagonal system may not be diagonally dominant. In this case 

3 

the inversion algorithm may lead to large error growth. Diag- 
onal dominance can be achieved by a spline adaptation of the 
finite-difference procedure given in Ref. 14. For all the 
problems treated here this modification is unnecessary. In other 
applications it will play a significant role if accurate solu- 
tions are to be obtained. 


VI, RESULTS 

* 

♦ 

Several model problems have been considered in order to 

evaluate the cubic spline collocation methods presented herein. 

For each of these problems an analytic solution or reliable 

numerical solution is available for comparison purposes. Spline 

interpolation (spline 2 and spline 4) is used to approximate the 

spatial gradients. For the one-dimensional Burgers equation the 

integration procedure outlined in Section II is adopted. Implicit 

* 

or two-step explicit methods are used. For the two-dimensional 
diffusion equation, solutions are obtained with the SADI formula- 
tion. The Laplace equation in cartesian and polar coordinates 
is evaluated with a spline successive approximation procedure. 
Finally, the similarity equations for the flat plate boundary-layer 
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and the two-dimensional stagnation point are solved by direct 
integration of the resulting two-point boundary value problems. 

Solutions are obtained with both uniform and non-uniform 
meshes. Three-point finite-difference calculations are included 
in order to assess the relative increase in accuracy associated 
with the higher-order procedures. The results are presented in 
tabular form so lhat meaningful comparisons arc possible. 

A. Burgers Equation 

The nonlinear Burgers equation (3a), with x^x, u=u(x,t) 
and x=x-(l/2)t, becomes 

w t - (u-l/2)u x =vu xx (8a) 

with v constant and the boundary conditions 

u-1 as x-» -co and iW) as x-*oo . (8b) 

The steady state solution of (8a) is 

u=[l-tanh(x/4v) 3/2 (8c) 

Spline 2 and the finite-difference solutions of (8) have 

“fi 

been discussed in Refs, 1 and 2. Both implicit and two-step 

explicit integration techniques, as outlined in Section IX, 

1 2 

have been applied successfully “ . Spline 4 solutions have now 
been obtained with the implicit and/or two-step procedures of 
(2a), (2b). The system (3) with the coefficients (6) are con- 

sidered. in the actual calculations the 3x3 system (3) is 

The nonlinear coefficient u is treated iteratively or with 
quasi-linearization 1 ' 2 . 


19 


reduced to a 2x2 systom. in is eliminated with (lo) or 

(If). The boundary conditions (Ob) on Uj, u^^ arc specified 

at x=x _ , with x >3, The boundary conditions on M. arc ob 
max max*~ i 

tainod from the third-order accurate relation 

(u xx , i-U“ (u x x ) i “ < 9i 


where i=l or N. 

The boundary condition (9a) can be applied , a two forms. 
These are outlined for the boundary i=l : 

(a) (u xx )l= ^ u xx^ 3 ” (M 3” Ml) 

With (u ) a evaluated from (5c) , we obtain 

(u xx )x=M 1 + (A/6) , (9b) 

♦ 

where a= h 3 /h 2 . From the governing equation (8a), 

* ♦ 

(u xx ) i = (u x ) 1 /2v= s m 1 /2 v 
so that with (If) , 

m 1 --h 3 M 1 /3~h s M JJ /6 + (u H -l)/h 0 


and (9b) becomes 

% ( 2 vt a v A/3+h s /3 ) -l-M 3 ( h 3 /6 - v ( 1+ c) A/3) t ( vA/3 ) M 3 -u g /h 3 1/h g . 

(9c) 

(b) An alternate form of (9a) , relating only the two points, 
i=l and 1=2, can bo derived by evaluating (u ) from (8a). The 

AX ^ 

temporal discretization is given by (2a) „ We obtain 


+ a„M" + a 3 u'„' = a, 


. M n+1 
33 


.n+l 


(9d) 
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whore 


! 


aj,« At (vh(u”*H0,5)h s /6) ? 
a a *=~At(v-(Ug-0<,25)h 3 /3) ; 

' 

a a = (l+At (u a -l)/h 3 ) ; 

a 4 «(u a + {u”-X) At/hJ . 

Per spline 4, A=(H*a a )/a(H‘cO 2 • For spline 2, sot A=0 so that 
(9c) is second-order accurate. Eq. (9b) is third-order accurate 
for both spline 2 and spline 4. Similar relations are obtained 
for the other boundary, where u N+1 f=0. 

ft 

The condition (9c) is independent of the time step At and 
somewhat less cumbersome. It was found that the accuracy of 
the solutions and the time to attain a converged steady state 
solution were virtually insensitive to the choice of the boundary 
condition (9c) or (9d) . This conclusion remains unchanged if 
the higher order effects in (9b), i.c., those terms multiplied 
by A, are treated explicitly in (9c) „ In this way (9c) reduces 
to a two-point implicit formula, in several cases the simpler 
spline 2 boundary conditions were applied with the spline 4 
procedure? the solutions always fell between the results of 
spline 2 and spline 4, but generally closer to those of spline 4. 
Therefore, if simplicity of boundary conditions is desired this 
is a reasonable approximation. 

Typical results are shown, for v=l/8, 1/16, 1/24 on Tables 
1-5. The increase in accuracy as one progresses from the finite- 
difference results to those of spline 2 and finally to spline 4 
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I 


is apparent. This is particularly true with the non-uniform 
meshes of Tables 2 and 3, For the conditions of Table 3, the 
finite-difference calculations with the two-step explicit pro- 
cedure did not converge, An oscillatory behavior was observed 
after 3200 iterations. In certain cases, where h^ is relatively 
large, the nature of the truncation errors (4a, 4b) of spline 2 
and spline 4 is such that a local value obtained with spline 2 
may be as accurate or more accurate than that obtained with 
spline 4. These are exceptional cases, however, and never occur 
for h.«l. A percentage error plot for the results of Table 1 
is shown on Figure 1. u e (x) denotes the exact solution (8c). 

Solutions for other v values are of a similar nature and 
therefore have not been included here. 

B. Linear Burgers Equation 

Consider the equation 

u x +vUx X =0 / on 0<x<l , 

with boundary conditions u(l)=l and on x=0, vu +u=0. The exact 

17 

solution is u e (x)=exp(l-x)/v. Fix has shown that with this 
flux boundary condition linear finite element theory naturally 
satisfies the required conservation condition at the boundary 
and therefore leads to more accurate solutions than obtained with 
non-divergence versions of spline 2 or conventional finite- 
difference theory o 

If finite-difference theory is developed in divergence or 
conservation form, the resulting equations are identical with 
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those of the linear, second-order accurate, finite element 

method. If spline 2 is recast in divergence form the solutions 

are considerably more accurate than the non- conservation results 

and also improve upon the conservation finite element (finite- 

difference) calculations. Therefore, with the flux boundary 

condition it appears that divergence form may be required if 

accurate spline solutions are to be obtained. On the other 

hand, if a modified derivative boundary condition was considered 

/ 

in lieu of the flux condition, the sensitivity to divergence form 
was no longer apparent. It is possible therefore that the flux 
condition represents a singular case. 

The governing systems of equations and the boundary condi- 
tions for the different formulations are as follows: 

Finite-Dif ference/Non-Divorgence Form 

v(ll i+j. +u i-r 2 u i >/h + (u i+r u i-i )/2 ° 0 (10a) 

At X=l, (10b) 

At x=0, v(u x - u )/2h + u =0 (10c) 

“ J- U 

Eliminating u_ x from (10c) with the difference equation (10a), 
we obtain 

2 v s (u 1 -u o )/h 2 H-u o (2v/h-l) = 0 . (lOd) 

Finite-Difference/Divergence Form = Finite Element 

(vu +u) = 0 

A 

Therefore, (vu^l-u) = ( v u x +u) 
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1 


or 


v < u iU +u i-l- 2u i )/h + < u i+l- u i-l>/ 2 “ 0 


(11a) 


The governing equation (11a) is identical with the non-divergence 
equation (10a) . The alterations appear in the boundary conditions 


At x=l , U N =1 

At -x~0, (vu +u). = 0, Therefore, 

X 

v(u 1 -u Q )/h + (uj+u^/2 = 0 . ’ (lib) 

The boundary condition (lib) differs frcft the non-divergence 

condition (lOd) 0 

Spline 2/Non -Divergence Form 

The governing equation (12a) is combined 

vM i +m i = 0 ' (12a) 

with the tpline relations (1). The boundary conditions are 
u N =l, vm o +u o = 0 * (12b) 

Spline 2/pivergence Form 

The governing equation (13a) is combined 

(vm+u) i+1 = (vm+u)^ (13a) 

with the spline relation (lc) « The boundary conditions are 

V x 

and 

vm +u +a( vmj+Ui ) -■= 0 , (13b) 

where a=0 corresponds to the exact boundary value and 

a=l corresponds to an averaged boundary condition. 
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The results of these calculations are shown on Table 6. 

It is seen that the non-conservation (NC) solutions with ten 

mesh points (N=10) are rather poor when compared with either 

the finite element or spline 2 conservation (c) solutions. It 

is significant that the spline divergence solutions, for both 

a=0 and a=l, are considerable improvements over the finite- 

element results o As the number of mash points increases the 

non-divergence solutions do show some improvement, with spline 2 

* 

more accurate than finite-differences, but these results are 
still less accurate than finite-element solutions. The ten point 
spline 2 divergence form ct=l solutions are about as accurate as 
the 50 point finite element results. 

Also shown on the table are ten point solutions with some- 
what modified derivative conditions at x=0. The exact solution 
is unchanged 0 ■ These derivative boundary conditions were treated 
in much the same manner as the flux condition for each of the 
procedures. For the finite-element solutions an average condition 
was applied. Significantly the large differences between diverg- 
ence solutions no longer occur . The spline solutions are always 
the most accurate, with a small increase in accuracy when diverg- 
ence form is assumed. 

C. Linear Corner Flow 

The two-dimensional diffusion equation 

u t = ^ (u xx +u yy > ' U=u(tffX ^ ) < 14a ) 

with the initial condition u(0,x,y)=0 and boundary conditions 
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u (t>0, 0,y>0) =1 , u (t>0, x>0, 0)-l 


u(t,x,y)->0 as x,y-*oo 


(14b) 


has the exact solution 


u=l~er£ X erf Y , (14c) 

where X = |(R e /t)^, Y = |(R e /t) i$ . 

This solution ’describes the impulsive motion of a right- 
angled corner formed by two infinite flat plates and results of 
the SADI spline 2 calculation have been presented in Refs. 1,2. 

The SADI procedure for the diffusion equation (14a) for both 
spline 2 and spline 4 is given as follows: 

Step 1, ujf - u" + ((u xx )^ +<u )" )At/(2R e ) US.) 


Step 2: +(t vv )?1' :L )At/(2R a ) 


IP 


XX' XJ 


yy i] 


where (« xx > L f M ij + <■ V 6) j" <1+a x> M ij +0 x M i-l. j ) 

and ( Vij = L ij +( V 6 )(L i.j+i-‘ (Wo y )L ij +a y L i-j-i ) 


(15b) 
. (16a) 
(16b) 


L. • and M. . each satisfy a tridiagonal equation of the form 

X J X J 

(lb) , 

A x * (l+oJ)/a x (ltcr x ) 2 ; A y = (l+aj)/a y (l+o y ) 2 


a x “ h i-H /h i ' °y “ k j+l A :i ; h i =x i“ x i-l ■ k j = yj-yj-l- 


The spline 2 formulation is recovered with A^A -0., The boundary 
conditions for u^ are given by (14b) . The boundary conditions 
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for L. .,M, . are obtained from (16) with (u )..»(u ) . ,«0 on 

13 id id yy 3-D 

the boundaries, or from (9a) with ( u xx )i + i # j obtained from (14a). 

The solution for step 1 is obtained with the tridiagonal 
2x2 system for M.. and u . . as described by (15a) and (lb). A 

X j 1 J 

similar procedure for L. *,u. . is required for step 2. 

1 J X J 

The solution for R q =1000 is given on Table 7. A non-uniform 
21x21 mesh with <^=0^=1. 5 was prescribed. The step size Ats=0.01. 
The solution is shown for t=2.0. All of the solutions are 

t 

reasonably good for this case, but once again the spline solutions 
are somewhat better. 

D, Laplace Equation 

The Laplace equation 


u +u - 0 7 u=u(x,y) , 
xx yy ' J 


(17a) 


with the boundary conditions u(0,y)=u(l,y)=0; u (x, 0) ^sinTTx? 
lim y-'oo u(x,y)=sO 

has the solution 

u(x,y) = (sinnx)exp(-ny) „ (17b) 

This boundary layer-like problem was chosen in order to 
evaluate the accuracy of the spline procedures, in particular 
with non-uniform meshes, when large gradients exist only over a 
limited region. In addition^ this problem will serve as a 
prototype for spline integration using successive approximation 
(relaxation) procedures. 


This procedure has been demonstrated for the Burgers equation 
by the discussion leading to (9d) . 
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Using the general expression for second derivatives ( 16 ) ^ 
Eq. (17a) can be put into a spline form* With the tridiagonal 
relationship for L ... and M . . (lb), this leads to a 3x3 system 

x j X j 

for the vector V. . at all interior mesh points: 

X j 

A i j v i, j-i +B i j v i j +c i j v i, 5+i +D i j +E ij v i+i, j = 0 < 18 > 


where 


ru . . 


V ij 


ij 


13 

M . . 
LiD 


A ij 


0 

-6/kj 


O A /6 

y y 

l 

o 


0 

0 

0 


B • * = 
13 


EJ 


0 ' 


(l^Oy) 


c 1+a 
6 x 

- - 

hr 


L i 


x 


l-i- a 

(1- -TT* 


1+c 

Ay) (1- “g A x ) 


2 (1+CTy) 


2 (1+a ) 
x' 


C ij 



A y /6 

0 


0 

* 

0 

a A 

X X 
6 

-6 


0 

7 D ij 

A 

0 ‘ 

0 

a 

y 3 

°y 

0 

0 

0 

0 


-6 

0 

1 
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E. ■ = 

3-D 


V 6 


-6 

o rnr 


The solution is obtained with a successive point relaxation 


procedure, 


n «ir A v oc+d , c v 0c) , D v (h+d _ v oo 1 

v ij - B ijL A ij v i, j-i h c ij v i, j+i 1 D ij v i-i-, j + E ij v i+i,j-i 


where the superscript k represents the iteration parameter. 

The system. is diagonally dominant and the eigenvalues of the 
amplification matrix, see Section V„ are all less than or equal 
to one. The results of this compuation are presented in Table 8, 
The values of u (.5,0) and u(.5,h_) are compared with the exact 

y 3 

solution (17b) . Also included in this table are the results 

obtained with the three-point finite-difference approximation 

for u and u . In order to make a more definitive comparison 
xx yy 

between the spline and finite-difference solutions, the surface 
value of u^ in the latter case was obtained by spline fitting 
the numerical values of u(x,y) . In one case noted on Table 8, 
a three-point end difference formula was applied. All of the 
calculations were performed with 10 mesh points in the normal 
or y-direction. In certain cases, spline 2 was used in the y- 
direction and spline 4 in the x-direction. These solutions are 
noted accordingly. 

The spline 4 results are the most accurate in all cases. 

For a uniform mesh the finite-difference and spline 2 results 
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are of equal accuracy as there are no convection effects in the 
problem. Moreover, if the spline 2 and the finite-difference 
solutions are averaged, the spline 4 results are closely approx- 
imated. For a non-uniform mesh the improved accuracy of spline 2 
over the finite-difference approximation is now apparent. 

The spline 4 results are remarkably accurate with o=1.7, 
h_=0.1 and y =28.66. For this mesh there are only four points 
in the region 0<y<l as compared with a uniform mesh (h=0.1) and 
ten points. The coarse mesh, spline 4 results are more accurate 
than the uniform mesh finite-difference solutions. 

i 

The 1.7/. 2 notation for a means that o=1.7 for Iv <0.2. 

For h.>0 o 2, a becomes unity. In this way the mesh width does 
not exceed a specified maximum value. This type of mesh alignment 
is useful in boundary layer problems, where a fine grid is desired 
near the surface, *nd a uniform but coarser mesh is required in 
the outer inviscid regions. This procedure is also applied for 
the boundary layer solutions in Section VI. F. An error plot is 
given on Figure 2 „ 

E. Potential Flov; Over a Circular Cylinder 

The governing equation in cylindrical coonrdinatcs for the 
potential flow over a circular cylinder is given by 


u + 
rr 




0 . 


( 20 ) 


The boundary conditions are u (1, 0)=O and lim u(r, 0) -rcosg. The 

r-too 

exact solution of Eq„ (20) with these boundary conditions is 
u=(r+-) cose, Eqso (17a) and (20) differ only by the appearance 
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of tho u term which is discretized toy the relations (le) or 

JU 

(If). The resulting 3x3 system for u.., L. . and M. . is of the 

Xj Xj 

form ( 19 ) , The coefficient matrices in tho present case will 
toe somewhat altered toy the u term. 

The results of the iterative solution are presented in 
Tables 9 - 11 . As in the previous examples, the finite-difference 
solutions are obtained toy using three-point central difference 
formulas. In Table 11, the slip velocity on the fore surface 
of the cylinder is presented. The superiority of the spline 
solutions over those resulting from finite-difference discretiza- 
tion is evident. It should toe noted that the slip velocity in 
the finite-difference case is obtained toy using a three-point 
central difference formula, while the spline solutions require 
only the two-point formula (le) . The higher accuracy of the 
two-point spline formula over the three-point finite-difference 
relations can be of considerable importance for problems witn 
derivative boundary conditions. 

F. Similarity Boundary Layers 

The boundary layer equations for the flow 
(3=0) and the two-dimensional stagnation point 
reduced to the following ordinary differential 
appropriate similarity transformations, 

u"H-fu / -(-3 (l-u 3 ) = 0 

f' = u 

The boundary conditions are 

f(o)=0, u(o)=0, lim u(x)=1.0 ' 

x-oo 


over a flat plate 
(3=1) can be 
system by using 

(21a) 

(21b) 

(21c) 
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Accurate numerical solutions have been reported in the litora- 

X 0 

ture [see Rosonhead xo ]. 

In the spline 2 and spline 4 formulation, Eq 0 (21a) is re- 


duced to a 2x2 system for u. and M. and the two-point boundary 

I X 


value problem is solved subject to (21c) . For the first-order 
equation (21b) , we obtain the following spline approximation 
from (If) : 


£ i+l* sf i +h i+l u i + 


h? . . 
H'l 


(MF^+ 


• 5 «W 


( 22 ) 


where MF. = (u')* for spline 2, For spline 4, the following 
relation to evaluate is easily derived from (le) and (If) : 

MP i+l +MF i= 2 (< u > U-l-M i)/ h i+l (23) 

MP 0 = (U ') 0 . 

Eqs, (22) and (23) give rise to an initial value problem for f ^ 
and MF^ w/iich is solved by a marching procedure. Eq. (22) 
leads to third-order accurate expression for therefore, for 

non-uniform meshes and third-order accurate solutions, this 
approximation is adequate even for spline -4. For the finite- 
difference solutions, a second-order accurate two-point formula 
for f^, which is consistent with the accuracy of the overall 
scheme, is obtained with the trapezoidal rule.' For (3=1, the 
nonlinear term u 2 is treated by quasilinearization so that 

(u *« )a = u 00p u (**>_. u oo, „ 

k is the iteration parameter. 
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The results of there computations for both uniform as 
well as non-uniform meshes are tabulated in Tables 12-16. The 
shear at the wall is proportional to f"(o) and this term has 
been evaluated for a variety of meshes. The results arc given 
on Tables 15 and 16 for (3=0 and 0=1, respectively „ Finite- 
difference solutions for u. are obtained by using the throe- 
point central difference approximation. 

As noted previously/, the notation o=l. 8/1 means that o=1.8 
until h. reaches 1.0, at which point c=1.0. h is tho first 
mesh width off the wall x=0; N is the total number of mesh 
points. N a is the number of mesh points in the boundary layer 
defined by x<6. At x=6, |u-1.0 |<10~ G , 

6=0 Blasius solution : 

The spline 4 solution for N=61, h a =0.1 and o=1.0 is almost 
identical with the "exact" solution of f " (o) =0 .469600 . ^ If 
spline 2 boundary conditions are used with a spline 4 interior 
point formulation, £"{o)=0. 46 9608. As previously noted, this value 
lies between the spline 2 and spline 4 results. With o=l . 8/2 , 
h 3 =0.5, N=21 and only 5 points in the boundary layer (N Q =5) , 
the spline 2 value of f"(o) is in error by only 2 %. For the 
larger h s values the spline 2 solutions are even more accurate 
than those found with spline 4„ Similar behavior was observed 
with Burgers equation in Section VI .A. An error plot is given 
on Figure 3 0 

6=1 stagnation point flow :- 

For 3=0, the exact solution has u'"(o) =u xv (o) =0 and there- 
fore the inherent lower-order accuracy of the finite-difference 
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calculation is somewhat obscured near the wall x«0, For the 
stagnation point solution where f "(o)«1.23258B, the improvement 
associated with the spline formulation is clearly demonstrated. 
Therefore, it would appear that spline integration should be 
extremely useful for boundary layor problems. 


G. Non-Similar Boundary Layer Analysis 

As a final test of the spline procedures the constant 
pressure boundary layer equations written in physical variables 
(x,y) were considered: 


uu + vu 
x y 




u 


X 


+ V. 


0 


The boundary conditions are 

y—0 : u = v = 0 
y»l u -♦ 1 


The initial conditions were given by 

u(o„y) = 1.0, y^O and u(o,o)=0 
v(o,y) = 0 

The equations were integrated for a Reynolds number Re=10 r ' and 
a non-uniform mesh of ten points normal to the surface. The 
solution for the normalized shin friction is shown on Figure 4. 
The value N 0 denotes the actual number of points within the 
boundary layer. The same criteria of Section VI. F was applied. 
As the boundary layor grows with distance x, N B increases. 

With 6 to 7 points in the final boundary layer profiles, the 
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spline 4 solutions aro quite accurate 


VI SUMMARY 

It has been demonstrated that higher-order calculation 
procedures using cubic spline collocation provide accurate 
solutions to a number of model problems,. The spline methods 
termed spline 2 and spline 4 can be used for two-point boundary 
value problems, as well as implicit, explicit, two-step, ADI 
and iterative integration procedures „ 

Spline 4 is fourth-order accurate with a uniform mesh and 
third-order with a moderate non-uniform mesh„ fc'pline 2 is 
second-order accurate for diffusion terms and fourth-order 
(third-order) for convection with a uniform (non-uniform) mesh. 
Derivative boundary values are obtained directly without the 
need for end differencing. For implicit linear systems, the 
spline methods remain unconditionally stable. 

The results confirm the higher-order accuracy of the spline 
methods and lead to the hopeful conclusion that accurate solu- 
tions for more practical flew problems can be obtained with 
relatively coarse non-uniform meshes. The relationship between 
the spline methods and similar Hermitian finite-difference pro- 
cedures must be explored. 

There has been no attempt to optimize the temporal inte- 
gration procedure so as to minimize computer times or increase 
temporal accuracy. The finite-difference calculations run 20% 
to 25% faster than the spline integrations. When spline fitting 
is used to evaluate finite-difference derivatives, as in 
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Section VI. c, the computer times are comparable. It is antici- 
pated that the reduced mesh requirements with these spline 
methods will result in a net improvement in computer storage 
and time. 
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TABLE I 


SOLUTION OF BURGERS EQUATION 
Z/ = l/8 J 1 (T = I.O, 31 EQUALLY 
SPACED POINTS 


§^9 

FD. 

SPLINE 2 

SPLINE 4 

EXACT 

0 

0.5000 

0.5000 

0.5000 

0.5000 

-0.2 

0.6999 

0.6860 

0.6900 

0-6900 

-0.4 

0.3447 

0.8290 

0.8322 

0.8320 

-0.6 

0.9269 

0.9IS0 

0.9170 

0.9 S70 

- 0.8 

0.9673 

0.9620 

0.9609 

0.9610 

- 1.0 

0.9857 

0.9830 

0.S820 

0.9820 

- 1.2 


0.9930 

0.9913 

0.9920 

- j,4 

0.9973 

0.9970 

0.9963 

0.9960 

- 1.6 

0.9988 

0.9990 

0.9983 

0.9980 

— 8.8 

0.9995 

0.9990 

0.9993 

0.3990 
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TABLE 2 

SOLUTION OF BURGERS EQUATION 
1 / = l/8,cr = 1.2, 15 POINTS 


gm 

ESsS 

F 0. 

SPLINE 2 

SPLINE 4 

EXACT 

0 

0.5000 

0.5000 

0.5000 

0.5000 

-0.3859 

0,9510 

0.8214 

0.8297 

0.8240 

-0.8^94 

o 

o 

o 

0.9778 

0.9654 

0.9676 

-6 .4080 

0.9990 

t .0004 

0.995 1 

0.9964 

-2.0750 

s.o 

1.0 

0.9989 

0.9S97 

-2.8770 

1.0 

6.0 

0.9999 

S.O 

-3.0420 

6.0 

6.0 

0.S99S 

6.0 

-5,000 

6.0 

6.0 

6.0 

6.0 
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TABLE 3 


SOLUTION OF BURGERS 
EQUATION! Z/=!/S,cr =1.8, 
15 POINTS 


u 

X 

SPLINE 2 

SPLINE 4 

EXACT 

0 

0.5000 

0.5000 

0.5000 

-0.0662 

0.5740 

mi 

0.5659 

-0.6855 

0.6S86 

0.6858 

0.6774 

-0.4001 

0.0695 

0.8452 

0.8321 

-0.7864 

6.0062 

0.9689 

0.9587 

-6 .4868 

1 .0665 

S .0083 

0.997*3 

-2.7354 

S .02ST 

1.0257 

6.0 

-4.9864 

6.0 

6.0 

8.0 
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TABLE 4 


SOLUTION OF BURGERS 

equation: v=\A6, <r=i.o, 

19 EQUALLY SPACED POINTS 


u 

X \ 

HD. 

SPLINE 2 

SPLINE 4 

EXACT 

0 

0.5000 

0.5000 

0.5000 

0.5000 

~ 0.2 

0.9000 

0.823 1! 

0.8356 

0.8320 

0. 4 

0.9878 

0.9046 

0.96 S 7 

0.9608 

- 0.6 

0.9980 

0.9952 

0.99 S 8 

0.99 SB 

0. © 

0.9998 

0.S995 

0.9982 

0. 9983 

- s.o 

1.0 

0.9999 

0.9996 

0.9997 

- S .2 

1.0 

S.O 

0.9999 

0.9999 
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TABLE 5 


SOLUTION OF BURGERS 
EQUATION! v = 1/24, <r =1.2, 
31 POINTS 


u 

x 

SPLINE 2 

SPLINE 4 

EXACT 

0 

0.5000 

0.5000 

0.5000 

-0.0688 

0.6936 

0.6957 

0.6955 

-0.8584 

0.86 68 

0.8S06 

0.8602 

-0.2505 

0.9506 

0.S526 

0.9529 

-0.3695 

0.9928 

0.9876 

0.9883 

-0.5622 

0.9995 

0.9975 

0.9979 

-0.6835 

6.0 

0.9990 

0.9997 

-0.S890 

6.0 

6.0 

6.0 

-6.6356 

6.0 

6.0 

6.0 

-4.9582 

6.0 

6.0 

6.0 
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EXACT SOLUTiON U Q = 148.4 , v = i/5 


FABLE 7 LINEARIZED CORNER FLOW 
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TABLE 8 

SOLUTION OF THE LAPLACE 
EQUATION 


METHOD 

-UJ.5,0) 

U 

c r 

ha 

^2 


EXACT SOLUTION 
FINITE DIFFERENCE 
SPLINE 2 

SPLINE 4X, SPLINE 2Y 
SPLINE 4 

3.142 

3.123 

3.164 

3.S54 

3.142 

0.7304 

0.7322 

0.7286 

0.7295 

0.7304 

1.0 

0.1 

0.1 

1.0 

EXACT SOLUTION 
SPLINE 4 

SPLINE 4X» SPLINE 2Y 
FINITE DIFFERENCE 

3.142 
3. 162 
3. 193 
2.82S* 

0.5335 

0.5329 

0.5277 

0.5401 

1.0 


0.2 

2.0 

EXACT SOLUTION 
FINITE DIFFERENCE 
SPLINE 2 

SPLINE 4X, SPLINE 2Y 
SPLINE 4 

3 . 142 
3. 189 
3. 175 
3 . S62 
3. 137 

0.9691 

0.9686 

0.9607 

0.9689 

0.9691 

1.7 


0.01 


EXACT SOLUTION 
FINITE DIFFERENCE 
SPLINE 2 

SPLINE 4X, SPLINE 2Y 

fcMrgjl 

El” 

0.939 1 
0.938 E 
0.9384 
0.9387 

S .7 

0.1 

0.02 

2.73 

EXACT SOLUTION 
FINITE DIFFERENCE 
SPLINE 4X, SPLINE 2Y 
SPLINE 4 

3 . 142 
3.220 
3. §05 
3 . 630 

0.7304 
0.7233 
0.7 26 3 
0.7313 

1.7 

I 



EXACT SOLUTION 
FINITE DIFFERENCE 
SPLINE 2 

SPLINE 4X, SPLINE 2Y 
SPLINE 4 

3.1416 
3.1694 
3.6676 
3J55S 
3. 1404 

0.2691 

0.9683 

0.9688 

0.9689 

0..9G96 

% 

OJ 

o.os 

LOSS 


* EVALUATED BY 3 - POINT END-DIFFERENCE FORMULA 


45 































TABLE 9 POTENTIAL FLOW OVER A CIRCULAR 
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TABLE II 


SLIP VELOCITY ON 
THE FRONT OF A 
CIRCULAR CYL 
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TABLE 15 f"(0) FOR BLASIUS EQUATION 
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FIG. i NON -UNEAR BURGERS EQUATION: * = l/8 



LAPLACE EQUATION 
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